\(\int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx\) [785]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [C] (warning: unable to verify)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [A] (verification not implemented)
   Giac [A] (verification not implemented)
   Mupad [F(-1)]

Optimal result

Integrand size = 31, antiderivative size = 120 \[ \int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2 a A x^{3/2} \sqrt {a^2+2 a b x+b^2 x^2}}{3 (a+b x)}+\frac {2 (A b+a B) x^{5/2} \sqrt {a^2+2 a b x+b^2 x^2}}{5 (a+b x)}+\frac {2 b B x^{7/2} \sqrt {a^2+2 a b x+b^2 x^2}}{7 (a+b x)} \]

[Out]

2/3*a*A*x^(3/2)*((b*x+a)^2)^(1/2)/(b*x+a)+2/5*(A*b+B*a)*x^(5/2)*((b*x+a)^2)^(1/2)/(b*x+a)+2/7*b*B*x^(7/2)*((b*
x+a)^2)^(1/2)/(b*x+a)

Rubi [A] (verified)

Time = 0.03 (sec) , antiderivative size = 120, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 2, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.065, Rules used = {784, 77} \[ \int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2 x^{5/2} \sqrt {a^2+2 a b x+b^2 x^2} (a B+A b)}{5 (a+b x)}+\frac {2 a A x^{3/2} \sqrt {a^2+2 a b x+b^2 x^2}}{3 (a+b x)}+\frac {2 b B x^{7/2} \sqrt {a^2+2 a b x+b^2 x^2}}{7 (a+b x)} \]

[In]

Int[Sqrt[x]*(A + B*x)*Sqrt[a^2 + 2*a*b*x + b^2*x^2],x]

[Out]

(2*a*A*x^(3/2)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(3*(a + b*x)) + (2*(A*b + a*B)*x^(5/2)*Sqrt[a^2 + 2*a*b*x + b^2*
x^2])/(5*(a + b*x)) + (2*b*B*x^(7/2)*Sqrt[a^2 + 2*a*b*x + b^2*x^2])/(7*(a + b*x))

Rule 77

Int[((d_.)*(x_))^(n_.)*((a_) + (b_.)*(x_))*((e_) + (f_.)*(x_))^(p_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*
x)*(d*x)^n*(e + f*x)^p, x], x] /; FreeQ[{a, b, d, e, f, n}, x] && IGtQ[p, 0] && (NeQ[n, -1] || EqQ[p, 1]) && N
eQ[b*e + a*f, 0] && ( !IntegerQ[n] || LtQ[9*p + 5*n, 0] || GeQ[n + p + 1, 0] || (GeQ[n + p + 2, 0] && Rational
Q[a, b, d, e, f])) && (NeQ[n + p + 3, 0] || EqQ[p, 1])

Rule 784

Int[((d_.) + (e_.)*(x_))^(m_.)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Dis
t[(a + b*x + c*x^2)^FracPart[p]/(c^IntPart[p]*(b/2 + c*x)^(2*FracPart[p])), Int[(d + e*x)^m*(f + g*x)*(b/2 + c
*x)^(2*p), x], x] /; FreeQ[{a, b, c, d, e, f, g, m}, x] && EqQ[b^2 - 4*a*c, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\sqrt {a^2+2 a b x+b^2 x^2} \int \sqrt {x} \left (a b+b^2 x\right ) (A+B x) \, dx}{a b+b^2 x} \\ & = \frac {\sqrt {a^2+2 a b x+b^2 x^2} \int \left (a A b \sqrt {x}+b (A b+a B) x^{3/2}+b^2 B x^{5/2}\right ) \, dx}{a b+b^2 x} \\ & = \frac {2 a A x^{3/2} \sqrt {a^2+2 a b x+b^2 x^2}}{3 (a+b x)}+\frac {2 (A b+a B) x^{5/2} \sqrt {a^2+2 a b x+b^2 x^2}}{5 (a+b x)}+\frac {2 b B x^{7/2} \sqrt {a^2+2 a b x+b^2 x^2}}{7 (a+b x)} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.01 (sec) , antiderivative size = 51, normalized size of antiderivative = 0.42 \[ \int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2 x^{3/2} \sqrt {(a+b x)^2} (7 a (5 A+3 B x)+3 b x (7 A+5 B x))}{105 (a+b x)} \]

[In]

Integrate[Sqrt[x]*(A + B*x)*Sqrt[a^2 + 2*a*b*x + b^2*x^2],x]

[Out]

(2*x^(3/2)*Sqrt[(a + b*x)^2]*(7*a*(5*A + 3*B*x) + 3*b*x*(7*A + 5*B*x)))/(105*(a + b*x))

Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 2.

Time = 0.08 (sec) , antiderivative size = 34, normalized size of antiderivative = 0.28

method result size
default \(\frac {2 \,\operatorname {csgn}\left (b x +a \right ) x^{\frac {3}{2}} \left (15 B b \,x^{2}+21 A b x +21 a B x +35 a A \right )}{105}\) \(34\)
gosper \(\frac {2 x^{\frac {3}{2}} \left (15 B b \,x^{2}+21 A b x +21 a B x +35 a A \right ) \sqrt {\left (b x +a \right )^{2}}}{105 \left (b x +a \right )}\) \(44\)
risch \(\frac {2 x^{\frac {3}{2}} \left (15 B b \,x^{2}+21 A b x +21 a B x +35 a A \right ) \sqrt {\left (b x +a \right )^{2}}}{105 \left (b x +a \right )}\) \(44\)

[In]

int((B*x+A)*x^(1/2)*((b*x+a)^2)^(1/2),x,method=_RETURNVERBOSE)

[Out]

2/105*csgn(b*x+a)*x^(3/2)*(15*B*b*x^2+21*A*b*x+21*B*a*x+35*A*a)

Fricas [A] (verification not implemented)

none

Time = 0.33 (sec) , antiderivative size = 30, normalized size of antiderivative = 0.25 \[ \int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2}{105} \, {\left (15 \, B b x^{3} + 35 \, A a x + 21 \, {\left (B a + A b\right )} x^{2}\right )} \sqrt {x} \]

[In]

integrate((B*x+A)*x^(1/2)*((b*x+a)^2)^(1/2),x, algorithm="fricas")

[Out]

2/105*(15*B*b*x^3 + 35*A*a*x + 21*(B*a + A*b)*x^2)*sqrt(x)

Sympy [F]

\[ \int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\int \sqrt {x} \left (A + B x\right ) \sqrt {\left (a + b x\right )^{2}}\, dx \]

[In]

integrate((B*x+A)*x**(1/2)*((b*x+a)**2)**(1/2),x)

[Out]

Integral(sqrt(x)*(A + B*x)*sqrt((a + b*x)**2), x)

Maxima [A] (verification not implemented)

none

Time = 0.21 (sec) , antiderivative size = 35, normalized size of antiderivative = 0.29 \[ \int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2}{35} \, {\left (5 \, b x^{2} + 7 \, a x\right )} B x^{\frac {3}{2}} + \frac {2}{15} \, {\left (3 \, b x^{2} + 5 \, a x\right )} A \sqrt {x} \]

[In]

integrate((B*x+A)*x^(1/2)*((b*x+a)^2)^(1/2),x, algorithm="maxima")

[Out]

2/35*(5*b*x^2 + 7*a*x)*B*x^(3/2) + 2/15*(3*b*x^2 + 5*a*x)*A*sqrt(x)

Giac [A] (verification not implemented)

none

Time = 0.27 (sec) , antiderivative size = 53, normalized size of antiderivative = 0.44 \[ \int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\frac {2}{7} \, B b x^{\frac {7}{2}} \mathrm {sgn}\left (b x + a\right ) + \frac {2}{5} \, B a x^{\frac {5}{2}} \mathrm {sgn}\left (b x + a\right ) + \frac {2}{5} \, A b x^{\frac {5}{2}} \mathrm {sgn}\left (b x + a\right ) + \frac {2}{3} \, A a x^{\frac {3}{2}} \mathrm {sgn}\left (b x + a\right ) \]

[In]

integrate((B*x+A)*x^(1/2)*((b*x+a)^2)^(1/2),x, algorithm="giac")

[Out]

2/7*B*b*x^(7/2)*sgn(b*x + a) + 2/5*B*a*x^(5/2)*sgn(b*x + a) + 2/5*A*b*x^(5/2)*sgn(b*x + a) + 2/3*A*a*x^(3/2)*s
gn(b*x + a)

Mupad [F(-1)]

Timed out. \[ \int \sqrt {x} (A+B x) \sqrt {a^2+2 a b x+b^2 x^2} \, dx=\int \sqrt {x}\,\sqrt {{\left (a+b\,x\right )}^2}\,\left (A+B\,x\right ) \,d x \]

[In]

int(x^(1/2)*((a + b*x)^2)^(1/2)*(A + B*x),x)

[Out]

int(x^(1/2)*((a + b*x)^2)^(1/2)*(A + B*x), x)